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ABSTRACT 

Aims. We study the Oxygen and Nitrogen abundances in the interstellar medium of high-redshift galaxies. 
Methods. We use high resolution and high signal-to-noise ratio spectra of Damped Lyman-a (DLA) systems detected 
along the hne-of-sight to quasars to derive robust abundance measurements from unsaturated metal absorption lines. 
Results. We present results for a sample of 16 high-redshift DLAs and strong sub-DLAs (log N{}1 l) > 19.5, 
2.4 < 2abs <3.6) including 13 new measurements. We find that the Oxygen to Iron abundance ratio is pretty much 
constant with [0/Fe] ~ -|-0.32±0.10 for —2.5 < [O/H] < —1.0 with a small scatter around this value. The Oxygen 
abundance follows quite well the Silicon abundance within ~0.2 dex although the Silicon abundance could be slightly 
smaller for [O/H] < —2. The distribution of the [N/0] abundance ratio, measured from components that are detected 
in both species, is somehow double peaked: five systems have [N/0] > —1 and nine systems have [N/O] < —1.15. 
In the diagram [N/O] versus [O/H], a loose plateau is possibly present at [N/O] ^ —0.9 that is below the so-called 
primary plateau as seen in local metal-poor dwarf galaxies ([N/O] in the range —0.57 to —0.74). No system is seen 
above this primary plateau whereas the majority of the systems lie well below with a large scatter. All this suggests 
a picture in which DLAs undergo successive star-bursts. During such an episode, the [N/O] ratio decreases sharply 
because of the rapid release of Oxygen by massive stars whereas inbetween two bursts. Nitrogen is released by low and 
intermediate-mass stars with a delay and the [N/O] ratio increases. 
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1. Introduction 

The production of Nitrogen in stars is the subject of strong 
interest as it is difficult to explain consistently the Nitrogen 
abundances measured in different astrophysical environ- 
ments. New abundance measurements in the interstellar 
medium (ISM) of our Galaxy as well as in external galaxies 
are of prime importance in these discussions. This could 
help decide what the sites of the Nitrogen production are. 
It is generally believed that the main contributors are the 
long-lived intermediate mass stars that are progenitors of 
Asymptotic-Giant-Branch (AGB) stars. The contribution 
from massive stars is still uncertain. 

The main nucleosynthetic pathway for the production 
of Nitrogen is the CNO cycle which takes place in the stel- 
lar H-burning layers. Nitrogen is thought to have both pri- 
mary and secondary origins depending on whether the seed 
Carbon and Oxygen nuclei are produced by the star itself 
(primary) or are already present in the interstellar medium 
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from which the star forms in which case Carbon and 
Oxygen seeds are left-overs from previous generations of 
stars (secondary) . Secondary production can happen in the 
H-burning layers of all stars as the Carbon seed is present, 
by definition, in these layers. Primary Nitrogen is produced 
when Carbon, synthetized in the Helium-burning shell of 
the star, penetrates into the Hydrogen-burning upper shell 
where it is transformed in Nitrogen by the CNO cycle. 
This happens in intermediate mass stars (4 < M/Mq < 7) 
during the Asymptotic Giant Branch phase (Henry et al. 
2000). Large uncertainties affect the theoretical predic- 
tions of both the primary and secondary nucleosynthesis of 
Nitrogen in low and intermediate-mass stars and the pos- 
sible contribution of primary Nitrogen from massive stars 
(Woosley & Weaver 1995, Marigo 2001, Maeder & Meynet 
2002). 

It is believed that the determination of the Nitrogen 
abundance in objects with metallicities spread over a large 
range may help answering this question. Therefore, it is 
usual to derive metallicities in nearby low-metallicity emis- 
sion line galaxies (e.g. Nava et al. 2006; Izotov et al. 2006) 
or low-metallicity stars in our Galaxy (Spite et al. 2005). 
A complementary approach is to derive metallicities di- 
rectly at high redshift in Damped Lyman-a (DLA) sys- 
tems (Pettini et al. 2002, Prochaska et al. 2002, Centurion 
et al. 2003) that have metallicities typically in the range 
—2.5 < Z/Zq < —1. Due to the large Hi column densi- 
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ties (log A^(H i) > 10^" cm"^) and the conspicuous pres- 
ence of metals, DLAs are believed to arise in high-redshift 
galaxies or at least to be located close to regions where 
star-formation occurs. 

It is usual to discuss these issues using the diagram giv- 
ing the Nitrogen to Oxygen abundance ratio versus the 
Oxygen abundance. In the case of secondary production, 
the ratio of Nitrogen to Oxygen abundances increases with 
increasing Oxygen mctallicity; whereas for primary produc- 
tion, the ratio remains contant as Nitrogen tracks Oxygen. 
In H II regions of nearby galaxies a trend in the [N/0] 
abundance ratio is seen with a primary plateau at low 
Oxygen abundances and a secondary behavior for abun- 
dances [0/H] > -1 (van Zee et al. 1998; Izotov & Thuan 
1999). DLA systems have a very different behavior in this 
diagram with [N/a] (see below) measurements well below 
the primary plateau (Pcttini et al. 2002; Centurion ct al. 
2003). This is probably related to the star- formation history 
of these objects. 

Oxygen and Nitrogen abundances can be derived di- 
rectly from the N{0 i)/N(R i) and N{N i)/N{R i) column 
density ratios. Because of efficient charge exchange reac- 
tions that link the neutral species together, the ionization 
corrections are negligible for log 7V(H i) > 19.5 (Viegas 
1995). Accurate O i and N i column densities are difficult 
to derive however. The main reason is that the absorp- 
tion lines arc located in the Lyman-a forest and are often 
blended. In addition, and this is especially true for O i, 
absorption lines are easily saturated. For N i the situation 
is less dramatic because the Nitrogen abundance is smaller 
and N i has two triplets around A1134 and A1200 A. On 
the contrary, the main O I absorption feature at A1302 is a 
single line and is often saturated. Other weaker absorption 
lines arc found much further in the blue, therefore deep 
in the Lyman-a forest and in a wavelength range where 
good signal-to-noise ratio (SNR) is difficult to obtain. This 
is why, in previous studies, the Oxygen abundance has of- 
ten been replaced by that of another a-element such as 
Sulphur or Silicon (Pettini et al. 2002, Centurion et al. 2003, 
Henry & Prochaska 2007) . Using Silicon instead of Oxygen, 
Centurion ct al. (2003) claimed that 75% of DLA systems 
show a mean value [N/Si] = —0.87 with 0.17 dex dispersion, 
corresponding approximately to the primary plateau ob- 
served locally, and 25% are clustered around [N/Si] = —1.5 
with even less dispersion (0.05 dex), the transition between 
low and high [N/Si] values happening at [N/H] = —2.8 (note 
that the latter transition may occur at [N/H] = —3, see 
Molaro et al. 2003). 

Here we present the (N/0) vs (0/H) diagram using 
robust determinations of Oxygen and Nitrogen mctallicities 
derived from unsaturated absorption lines selected from the 
largest UVES sample of DLA systems available up to now 
(Ledoux et al. 2003, 2006a). In Section 2 we present our 
sample. We describe the O and N abundance measurements 
in Section 3 and discuss the results and implications in 
Sections 4 and 5. 

2. UVES DLA sample 

Most of the systems in our sample were selected from 
the follow-up of the Large Bright QSO Survey (Wolfe et 
al. 1995) and observed at the VLT with UVES between 

2000 and 2004 in the course of a systematic search for 
molecular hydrogen at Zabs > 1-7 (Petitjean et al. 2000; 



Ledoux et al. 2003). Our sample comprises 61 bona-fide 
DLA systems (logAf(Hi) > 20.3) and 13 strong sub-DLA 
systems with total neutral hydrogen column densities in 
the range 19.5 ^ logiV(Hi) < 20.3. Characteristics of the 
systems (mctallicities, N(}i i) column densities, kinemat- 
ics), except for four sub-DLAs with log N(R i) = 19.7- 
19.8), are given in Ledoux et al. (2006a). The absorp- 
tion line analysis was performed in an homogeneous man- 
ner using standard Voigt-profile fitting techniques adopting 
the oscillator strengths compiled by Morton (2003). Total 
column densities were derived as the sum of the column 
densities measured in individual components of the line 
profiles. Average gas phase metallicities relative to solar, 
[X/H] = log[iV(X)/Ar(H)] -log[N{X)/N{R)]Q, were calcu- 
lated using solar abundances from Lodders (2003). 

As we are interested in measuring the [0/H] and [N/H] 
abundances from the A^(0 i), A^(N i) and A^(H i) column 
densities, we have to select systems where we can derive 
accurate column densities, e.g. where at least one transi- 
tion from both N I and O I is not strongly saturated. Note 
that the transitions from these elements are most of the 
time redshifted in the Lyman- a forest and we require in 
addition that at least one unsaturated transition be free 
of blending. From the above sample and under these condi- 
tions, we selected 13 sytems where we could measure A'^(0 i) 
and A''(N i) column densities accurately. 

Prom the literature, we added three measurements at 
^a^bs = 3.390, 2.844 and 3 toward, respectively, Q 0000-263 
(Molaro et al. 2001), Q 1946+769 and QXO 0001 
(Prochaska et al. 2002). As a consequence of our selection 
criteria, we did not include the system at 2abs = 

4.466 to- 
ward Q 0307-4945 (Dessauges-Zavadsky et al. 2001) as the 

complex O i and N i velocity profiles, spread over more than 
200 km s^^, are highly saturated and/or blended. Neither 
did we include the systems at Zabs = 2.076 and 2.456 to- 
ward, respectively, Q 2206-199 and Q 1409-J-095 (Pettini 
et al. 2002) because the Oxygen abundance can only be ill- 
defined from strongly saturated O iA1302 absorption lines. 
In the case of Q 2206—199, the strongest N i feature is also 
blended (Molaro et al. 2003). 

Ionization may be a concern when deriving the (N/0) 
abundance ratio from the N(N l)/N(0 l) column density 
ratio (Viegas 1995, Prochaska et al. 2002). If the gas is 
neutral then charge exchange reactions are fast enough so 
that O I and N l are both tied to H I and there is no need 
for correction for log iV(H i) > 19.5 (Viegas 1995). In case 
of the presence of enough hard photons (or cosmic rays) 
the ionization balance could be displaced towards higher 
ionization species and both O l and N l could be over- 
ionized compared to H i. However for most DLA systems, 
hard photons come predominantly from the background 
ionization field. The ionization parameter for this field at 
z = 3 is U < 10~^ for typical densities expected for DLAs 
(n > 1 cm~^; see e.g. Petitjean et al. 1992), and ioniza- 
tion corrections are again small. Indeed, Prochaska et al. 
(2002) inspected about twenty DLAs and concluded that 
the ionization correction for [N/a] is at most of the order 
of -1-0.1 dex. This amount should probably be added to the 
measurement errors. In the following we did not apply any 
correction as we show that there is no correlation between 
[N/0] and log iV(H l). 
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Name 


n 


^ahs 


log iV(H I) 

(cm 2) 


log N(0 i) 

(cm 2) 


cr 


[O/HJ 


log A^(N i) 
(cm 2) 


a 




"TnTTttT 

[N/HJ 


Q 0102-190 


1 


2.92646 


20.00+0.10 


14.68 


0.16 




12.67 


0.11 








2 


2.92661 




14.66 


0.17 




<12.35 


d 








3 


2.92727 




14.20 


0.02 




12.38 


0.09 








4 


2.92771 




14.41 


0.02 




12.55 


0.06 








Tot 






15.13 


0.08 


-1.56 


12.76 


0.11 




-3.07 


Q 0112-306 


1 


2.41844 


20.50+0.08 


14.78 


0.12 




12.66 


0.08 








2 


2.41861 




14.45 


0.06 




12.99 


0.04 








Tot 






14.95 


0.08 


-2.24 


13.16 


0.04 




-3.17 


Q 0347-383" 


1 


3.02463 


20.73+0.05 


16.12 


0.12 




14.15 


0.03 








2 


3.02485 




16.18 


0.18 




14.47 


0.03 








Tot 






16.45 


0.11 


—0.97 


14.64 


0.02 




— 1.92 


(4 0841+129 


1 


2.47622 


on on_i_n 1 n 

20.80+0.10 


16.09 


0.03 


— 1.40 


13.94 


0.02 




—2.69 


Q 0913+072*" 


1 


2.61829 


20.35+0.10 


14.08 


0.02 




<12.10 


d 








2 


2.61844 




14.43 


0.02 




<12.20 










Tot 






14.59 


0.02 


—2.45 


<12.45 




< 


—3.73 


Q 1108-077 


1 


3.60767 


20.37+0.07 


15.37 


0.03 


-1.69 


<12.84 


d 


< 


-3.36 


Q 1337+113 


1 


2.50793 


20.12+0.05 


14.90 


0.10 


-1.91 


12.84 


0.10 




-3.11 


Q 1337+113 


1 


2.79557 


21.00+0.08 


15.21 


0.12 




<12.40 










2 


2.79581 




15.59 


0.08 




13.98 


0.03 








Tot 






15.74 


0.07 


-1.95 


13.99 


0.03 




-2.84 


Q 1340-136 


1 


3.11835 


20.05±0.08 


15.52 


0.02 


-1.22 


13.28 


0.02 




-2.60 


Q 1409+095= 


1 


2.45595 


20.53+0.08 


15.05 


0.02 




<13.26 


d 








2 


2.45644 




15.00 


0.02 




<13.14 










Tot 






15.33 


0.02 


— 1.89 


<13.51 




< 


—2.85 


1 /inn i nnc^c 


1 


(i^^QO^ 
Z.DDozl 


1 n on_Ln no 


1 c 1 


n no 


— i.dl 


13. oU 


n no 




— 2.io 


Kci zuoy — oDU 


1 
1 




on QO-|_n ns 


10. Do 


U.UD 




AO 


n no 








9 


o.Uozyo 




iO.oO 


u.uo 




10. / 1 


u.uz 








3 


O.VJOiJ J-U 




15.02 


0.12 




12.85 


0.14 








Tot 






16.09 


0.04 


-1.58 


13.95 


0.02 




-2.86 


Q 2332-094 


1 


3.05633 


20.50+0.07 


14.97 


0.02 




<12.93 


d 








2 


3.05658 




14.46 


0.03 




<12.93 


d 








3 


3.05677 




14.71 


0.07 




<12.93 


d 








4 


3.05690 




15.32 


0.02 




<12.93 


d 








5 


3.05723 




15.67 


0.04 




13.73 


0.03 








6 


3.05737 




14.66 


0.04 




<12.93 


d 








Tot 






15.95 


0.02 


-1.24 


13.73 


0.03 




-2.60 



"also Ledoux et al. (2003); also Erni et al. (2006); " also Pettini et al. (2002); 3a detection limit. 



3. Column Densities and Abundances 

We performed in a usual manner Voigt profile fitting of 
absorption lines associated to the 13 systems we selected 
from our DLA sample (sec e.g. Ledoux ct al. 2006a, Erni et 
al. 2006). The overall decomposition in subcomponents is 
derived from the simultaneous fit of numerous absorption 
lines from different species (most importantly Zn II, Fe il 
and Si 11) redshifted outside the Lyman-a forest and free 
from blending. Wc used the resulting component structure 
to detect N i and O l features and confirm that they are not 
blended. Column densities obtained for each of the compo- 
nents arc given in Tabic 1. Fits arc shown in Figs. 8, 9 
and 10. The fits to the absorption lines in the Zabs — 3.025 
system towards Q 0347—383 can be seen in Ledoux et al. 
(2003). Errors for individual components are the 1 a errors 
from Voigt profile fitting as given by FitLyman. 
Following the usual procedure, [N/H] and [0/H] abun- 
dances are obtained by adding the column densities in all 
detected individual components. For O i, all the systems 
are dominated by one strong component in which we do 
not expect much ionization correction (see next Section for 
details on each of the systems). It is possible that in some 
of the weak components the gas is partially ionized and 



that some ionization correction should be applied to these 
components. However, this ionization correction is probably 
smaller than 0.2 dex in the weakest components and these 
components do not account, in total, for more than ~20% 
of the total column density. We therefore do believe that 
the error on [0/H] due to the procedure is less than 0.2 dex 
for the whole system (see also the discussion in Prochaska 
et al. 2002). For N i, the above error is less as the transi- 
tions are weaker and the species is always detected in less 
components. 

For the [N/0] ratio, we use the cohmm densities in the 
components where both species are detected. Note that 
when N I is detected in only one component, this is al- 
ways the strongest O i component. It is interesting that, 
in case N I and O I are detected in several components, 
as toward Q 2332—094, the [N/0] ratios are similar in all 
components. 

We adopt the Oxygen and Nitrogen solar abundances 
from Fodders (2003): 12 + (0/H)q = 8.69 and 
12 + (N/H)0 = 7.83. Note that the main discussion on 
the production of these elements is little dependent on the 
exact adopted values. Depletion onto dust-grains is known 
to be much smaller in DLA systems compared to the ISM 
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of our Galaxy mostly because of smaller (by a factor of 
at least ten) metallicities. Depletions should be therefore 
much smaller than in the diffuse galactic ISM clouds where 
it has been shown that, apart from local effects, the [N/0] 
ratio is similar to the solar ratio (Knauth et al. 2006). 

4. Comments on individual systems 

Q 0000—263: We adopt the column densities from Molaro 
et al. (2001). The Oxygen abundance is derived from unsat- 
urated O iAA925,950 absorption lines. The Nitrogen abun- 
dance is derived from the unsaturated A953 line and the 
A1134 triplet. The system is modeled with only one com- 
ponent. 

Q 0102—190 (Fig. 8, top panel): The structure of the system 
is derived from the fit of Fe ii and Si ii lines. A^(N i) is 
derived from a A1199 optically thin feature corresponding 
in redshift with the strongest O i component. We add all O i 
column densities (obtained from optically thin AA976,1039 
features) to calculate the O abundance. Note however that 
the red satellites do not contribute much. The [N/0] ratio 
is taken as the column density ratio in the main central 
component. 

Q 0112—306 (Fig. 8): The two-component structure of the 
system is very well defined from Si ii and Fc ii absorption 
lines. This is why we are confident that the log A^(0 i) col- 
umn density derived from the moderately saturated A1302 
absorption is robust. N i is detected through an optically 
thin All 99 feature. 

Q 0347—383: This system has been studied by Levshakov 
et al. (2002) and Prochaska et al. (2002). The UVES data 
have also been reanalysed by Ledoux et al. (2003). N(N i) is 
constrained by optically thin A1134 features and N{0 i) by 
A950 and A974. Note that Levshakov et al. (2002) derived 
log iV(N i) = 14.89, 0.25 dex larger than our result. We 
believe this is because the latter authors do not restrict 
their fit to optically thin lines. 

Q 0841-1-129 (Fig. 8): This is a beautiful one-component 
system with several optically thin N i components (see also 
Centurion et al. 2002). A^(0 i) is derived from the con- 
sistent fit of optically thick A1302 and 1039 features and 
the moderately saturated A950 absorption. Centurion et al. 
(2003) found a N i column density 0.16 dex larger. This is 
likely due to the difference in data quality. 
Q 0913-f072 (Fig. 8, bottom panel): The absorption pro- 
file is made of two closely blended components of similar 
strength. An upper limit on log N(N i) is obtained from 
A1199 (a feature is present but below the 3ct detection limit, 
see also Erni et al. 2006). Upper limits are derived from the 
noise in the adjacent continuum; when two components are 
present we conservatively consider each component inde- 
pendently. N{0 I) is derived from optically thin A1039. 
Q 1108—077 (Fig. 9, top panel): This is a one-component 
system. An upper limit on A^(N i) is obtained from A1134. 
A''(0 l) is derived from an optically thin A1039 feature. 
Q 1337-^113, Zabs = 2.508 (Fig. 9): This is again a single 
component system with well constrained N i and O i column 
densities from optically thin, respectively, A1200 and A1039 
absorption lines. 

Q 1337-^113, Zabs = 2.796 (Fig. 9): There are two compo- 
nents in this system. iV(N i) and N{0 i) are both well de- 
fined from optically thin, respectively, A1134,1200 and A976 
features. The component at Zabs = 2.79581 is the strongest 
and is the only component detected in N i. It is possible that 



N I is somewhat ionized in the component at Zabs = 2.79557. 
As done for all systems, we therefore add the O i column 
densities to calculate the Oxygen abundance but take the 
N/0 ratio from the component at z^bs — 2.79581. 
Q 1340—136 (Fig. 9, bottom panel): This is a beautiful one- 
component system with optically thin transitions. 
Q 1409-1-095, Zabs = 2.456 (Fig. 10, top panel): Although 
there may be some absorption at the wavelengths corre- 
sponding to the N I transitions, we derive only an upper 
limit on A^(N i) in this two-component system (see also 
Pcttini et al. 2002). The O i column density is well defined 
by optically thin A1039 absorption. 

Q 1409-1-095, Zabs = 2.668 (Fig. 10): The consistency of the 
AA1199,1200 and A1134 features makes us believe that N i 
is present in this one-component system. The O i column 
density is well defined by optically thin A1039 and A950 
absorptions. 

Q 1946-f 769: We use the measurements by Prochaska et al. 
(2002). 

Q 2059-360 (Fig. 10): N i and O i are detected in the 
three components of this system in, respectively, optically 
thin A1134,1200 and A950,976 transitions. It is important 
to note that the [N/0] ratio is similar for the three com- 
ponents. This supports the assumption made in this paper 
that ionization corrections are negligible. 
Q 2332—094 (Fig. 10, bottom panel): The system is com- 
plex but dominated by two main components. We add all 
O I column densities but the contribution of the three satel- 
lite components is negligible. N i is detected in a single 
component corresponding to the strongest O i component 
at Zabs — 3.05723 and we use the [N/0] ratio in this com- 
ponent. 

QXO 0001: We use the measurements by Prochaska et al. 
(2002). Note that only one digit is given by these authors 
for the absorption redshift. 

5. Results 

A summary of the abundances of interest here, measured 
in the sixteen systems, is given in Table 2. 

5.1. The Oxygen abundance in DLAs 

Few Oxygen abundance measurements have been reported 
up to now in DLAs (see Molaro et al. 2000; Pettini et al. 
2002). The main reason is that most of the time O iA1302 
absorption lines are badly blended or saturated. High reso- 
lution and high SNR spectra covering the bluer part of the 
optimal wavelength range are needed to detect weaker O I 
transitions. Thanks to the large number of DLAs observed 
in the course of the survey for molecular hydrogen and to 
the UVES sensitivity into the blue we could built a sample 
of robust measurements. 

Fig. [H gives the [0/Fe] abundance ratio versus the 
Oxygen abundance. In principle, this ratio should be cor- 
rected for depletion of metals onto dust-grains. However 
this correction is known to be negligible for [X/H] < —1.2 
(Prochaska & Wolfe 2003; Ledoux et al. 2003; Wolfe et al. 
2005). To support this idea, we plot in Fig.[2]the abundance 
ratio [S/Fe] versus the Sulfur abundance [S/H] for all sys- 
tems in our DLA sample (Ledoux et al. 2006a). Sulfur and 
Iron are known to be, respectively, little and strongly de- 
pleted onto dust-grains. It is apparent that the scatter of 
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Name 




[Fe/H] 


[O/H] 


[N/H] 


(0/H)+12= 


(N/O)" 


[S/H] 


[Si/H] 


Q 0000-263'^ 


3.390 


-2.00(0.09) 


-1.67(0.13) 


-2.50(0.08) 


7.02 


-1.69(0.10) 


-1.89(0.09) 


-1.88(0.08) 


Q 0102-190 


2.926 


-1.68(0.10) 


-1.56(0.13) 


-3.07(0.15) 


7.13 


-2.30(0.16) 


-1.37(0.10) 


-1.50(0.10) 


Q 0112-306 


2.418 


-2.64(0.09) 


-2.24(0.12) 


-3.17(0.09) 


6.45 


-1.79(0.09) 


-1.25(0.08) 


-2.42(0.08) 


Q 0347-383 


3.025 


-1.85(0.05) 


-0.97(0.12) 


-1.92(0.05) 


7.72 


-1.81(0.11) 


-1.18(0.06) 


-1.49(0.06) 


Q 0841+129 


2.476 


-1.79(0.10) 


-1.40(0.10) 


-2.69(0.10) 


7.29 


-2.15(0.04) 


-1.55(0.10) 


-1.58(0.10) 


Q 0913+072 


2.618 


-2.72(0.10) 


-2.45(0.10) 


< -3.73 


6.24 


< -2.14 


-1.66(0.10) 


-2.57(0.10) 


Q 1108-077 


3.608 


-1.96(0.07) 


-1.69(0.08) 


< -3.36 


7.00 


< -2.53 




-1.57(0.07) 


Q 1337+113 


2.508 


-2.23(0.05) 


-1.91(0.11) 


-3.11(0.11) 


6.78 


-2.06(0.14) 




-1.84(0.06) 


1 QQ7 1 1 1 Q 

ioor+iio 




— 2.i4^U.Uo ) 


— i.yo^U.iUj 


— 2.84^U.U9 ) 


0.(4: 


1 ^iT f r\ r\Ck\ 

— i.onu.uyj 




— i.o4(U.U9j 


Q 1340-136 


3.118 


-1.59(0.08) 


-1.22(0.08) 


-2.60(0.08) 


7.47 


-2.24(0.03) 


-1.41(0.08) 


-1.18(0.08) 


Q 1409+095 


2.456 


-2.43(0.08) 


-1.89(0.08) 


< -2.85 


6.80 


< -1.82 


-1.74(0.08) 


-2.17(0.08) 


Q 1409+095 


2.668 


-1.51(0.08) 


-1.31(0.08) 


-2.13(0.08) 


7.38 


-1.68(0.03) 


-1.40(0.09) 


-1.28(0.08) 


Q 1946+769'' 


2.844 


-2.50(0.06) 


-2.14(0.06) 


-3.51(0.07) 


6.55 


-2.23(0.04) 




-2.21(0.06) 


Q 2059-360 


3.083 


-1.97(0.08) 


-1.58(0.09) 


-2.86(0.08) 


7.11 


-2.14(0.04) 


-1.76(0.09) 


-1.66(0.09) 


Q 2332-094 


3.057 


-1.60(0.07) 


-1.24(0.07) 


-2.60(0.08) 


7.45 


-1.94(0.05) 


-1.32(0.08) 




QXO 0001* 


3 


< -1.08 


-1.62(0.05) 


-3.22(0.06) 


7.07 


-2.46(0.05) 




-1.79(0.05) 



"Molaro et al. (2001); "Prochaska et al. (2002); V = cr([0/H]) 

''Errors are taken as square root of the geometrical mean of the variances of A''(N l) and A'^(0 l) as given in 
Table 1 except for Q 0102-190, Q 1337+113 at z^bs = 2.796 and Q 2332-094 for which we consider only 
the main components of the systems to calculate the [N/0] ratio (see Text in Section 3). 




-2.0 -1.0 
[O/H] 




-2.0 -1.0 
[S/H] 



Fig. 1. [0/Fe] abundance ratio versus Oxygen abundance Fig. 2. Sulfur to Iron abundance ratio, [S/Fe], versus Sulfur 
[O/H]. The dashed and dotted lines are drawn for conve- abundance relative to solar, [S/H]. The lines are drawn for 



nience. 



convenience. 



the [S/Fe] values is much less below [S/H] = —1.2 than 
above. This behavior is a consequence of the largest deple- 
tion of Iron onto dust-grains for [S/H] > —1.2. 

It can be seen that because of our choice to reject 
the systems where all O i transitions are badly saturated, 
our sample is biased against metallicities [O/H] > — 1. 
All our systems except one have metallicities in the range 
-2.5 < [O/H] < -1.2. The stricking result here is that 
the [0/Fe] ratio is pretty constant over the range studied. 
The mean value of the 14 measurements at [O/H] < — 1 
(therefore excluding Q 0347-383) is [0/Fe] = 0.32+0.10 
with a slight tendency of this ratio to increase at lower 
metallicity. The only departure from this picture is the 
point at [O/H] = -0.97, [0/Fe] = 0.88, from the system 
at Zabs = 3.025 towards Q 0347—383. This is probably due 
to depletion onto dust-grains ([Zn/Fc] = 0.72). Actually, 



molecular hydrogen, although with very low molecular frac- 
tion, is detected in this system (Levshakov et al. 2002, 
Ledoux et al. 2003), supporting the presence of dust. 

As the Oxygen abundance is difficult to measure, the 
Si abundance is often used instead. Prochaska & Wolfe 
(2002) and Dessauges-Zavadsky et al. (2006) have found 
that the [Si/Fe] ratio is pretty constant amongst DLA sys- 
tems around a value of +0.43 with a relatively small dis- 
persion. We plot in Fig. [3] the [O/H] abundance versus the 
[Si/H] abundance. It is apparent that the two abundances 
correlate well. There may be a slight tendency for Oxygen 
to be more abundant than Silicon for [O/H] < —2 but there 
are only four measurements there. 
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Fig. 3. Oxygen abundance relative to solar, [0/H], versus Fig. 4. [N/0] abundance ratio versus neutral hydrogen col- 
Silicon abundance [Si/H]. Dotted line is for [0/H] = [Si/H]. unin density log N{}i i). The vertical dotted line indicates 

the standard definition of DLA systems (log iV(H i) > 20.3). 



5.2. The [N/0] abundance ratio 

Since we will use two different notations, we would 
like first to clarify them. The (N/0) abundance ratio 
is simply the ratio of the Nitrogen to Oxygen abun- 
dances, (0/H) = log7V(0 i)-logA^(H i) and the same 
for (N/H). This ratio is often plotted versus the Oxygen 
abundance written as (0/H) 4- 12. It is however very 
convenient, and usual, to refer to solar abundances, 
[0/H] = (0/H) — (O/H)©. In that case solar abundances 
should be specified. We use: 12 -I- (O/H)© = 8.69 and 
12 + (N/H)© = 7.83. Therefore [N/0] = (N/0) + 0.86. 
What is most important in the following is the value 
assigned to the so-called "primary plateau". Nava et 
al. (2006) observationally derive a primary plateau at 
(N/0)pp — —1.43 with objects distributed within a range of 
-1.54 to -1.27. In that case, [N/0]pp = -0.57. From mea- 
surements in blue compact dwarf (BCD) galaxies Izotov & 
Thuan (2004) however find slightly smaller values with a 
plateau at about (N/0)pp = -1.6 or [N/0]pp = -0.74. 

In Fig. d] the [N/0] abundance ratio is plotted versus 
the neutral hydrogen column density. It can be seen that 
there is no apparent correlation between the two parame- 
ters. The smallest [N/0] value is found for a system with 
log iV(H i) > 20.3 and one of the largest values is found 
for the system with the lowest log A''(H i). This gives ad- 
ditional confidence in the assumption that the ionization 
correction is negligible when deriving the abundances from 
A^(0 i), A^(N i) and iV(H i) (Viegas 1995). 

There is a possible dichotomy between systems with 
[N/0] -0.9 and -1.4. It can be seen on Fig. [5] that the 
distribution of the [N/0] abundance ratio has two peaks, 
one around [N/0] ^ —1.35 and one around [N/0] ^ —0.9. 
This could be however a consequence of small number 
statistics. Out of 13 systems with definite measurements, 
7 have [N/0] < -1.15 and 5 have [N/0] > -1.0. If we 
add the two upper limits with [N/0] < —1.2, we find that 
the number of systems with [N/0] < —1.15 is about twice 
larger than those with [N/0] > —1.0. A similar dichotomy 
has already been noted in the [N/a] ratio by Centurion et 
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in 



B 




o 
6 
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Q I \ I : \ I ] : ] I \ : I ] I \ ] 

-1.5 -1.0 
[N/0] 

Fig. 5. Histogram of [N/0] abundance ratio for the 13 sys- 
tems with definite measurements in our sample. 

al. (2003). They find however that 75% of systems clus- 
ter around [N/a] ^ —0.87 and the remaining 25% clus- 
ter around [N/a] ~ —1.45. Also, note that the dichotomy 
present in our data is less pronounced than what is claimed 
by Centurion et al. (2003) because of a large dispersion of 
the values below [N/0] ~ —1.2 (see also Henry & Prochaska 
2007) . One could wonder if these differences could be due 
to our sample being biased against high Oxygen metallic- 
ity systems compared to the Centurion et al. sample. We 
do not think this is the case because there is no obvious 
correlation between [0/H] and [N/0]. 

In addition. Centurion et al. (2003) claimed that there 
is a transition between the two regimes at [N/H] ~ —2.8 
with most, if not all, of the [N/a] < —1.5 systems hav- 
ing [N/H] < -2.8. Molaro et al. (2003) place this transi- 
tion at [N/H] ^ —3.0 and argue that this could be under- 
stood as the point at which Nitrogen production from AGB 
stars begins to dominate that of massive stars. Although 
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Fig. 6. Nitrogen to Oxygen abundance ratio, [N/0], versus 
Nitrogen abundance [N/H]. [X/H] = (X/H) - (X/H)©. 



all the points with [N/0] < -1.2 have [N/H] < -2.4, we 

find two systems with [N/0] > -1.2 and [N/H] 3.0. 

Again, below [N/0] < —1, the scatter in the measurements 
is larger than above this limit. Note however that the four 
points with [N/H] < -3.2 have [N/0] < -1.4 so that the 
figure is not inconsistent with some kind of break below 
[N/H] = —3.2. More data are needed below this limit be- 
fore any firm conclusion can be drawn. 

The classical plot giving [N/0] versus [0/H] is shown 
in Fig. [71 Several features have to be noted here. First, all 
the DLA measurements are located in the region delineated 
by the usual primary and secondary lines. Recall that, by 
definition, if Nitrogen is a primary element, it should be 
produced at the same time as Oxygen and the [N/0] ratio 
should be a constant whatever the Oxygen metallicity is. 
The primary plateau observed locally is derived from mea- 
surements of abundances in nearby low-metallicity galaxies: 
[N/0]pp in the range —0.57 to —0.74. If Nitrogen is a sec- 
ondary element, its production is favored at large Oxygen 
abundance and the [N/0] ratio should increase with in- 
creasing Oxygen abundance. The fact that all the DLA sys- 
tems are found in this region may indicate that at least part 
of the systems are in the transition zone as a consequence 
of delays in the release of heavy elements from intermediate 
mass stars. This is what would be expected if during a star- 
burst, at the beginning, high mass stars eject material with 
low [N/0] ratio and, after some time, intermediate mass 
stars eject material with higher ratio. The pathway of one 
particular system in the diagram during his lifetime could 
be a line starting from the left-bottom corner of the Figure 
towards the up-right direction. 

Secondly, DLA measurements are all below the local 
primary plateau ([N/0]pp ~ -0.57 to -0.74). If we be- 
lieve that the five top- most [N/0] DLA measurements de- 
fine what can be considered as a plateau, this plateau is at 
[N/0] ~ —0.9, therefore at least —0.15 dex below the local 
primary plateau. Note that this has not been recognized 
by studies based on [N/Si] or [N/S] ratios (Centurion et al. 
2002, Pettini et al. 2007). We understand this conclusion is 
based on small number statistics and should be confirmed 



with more data. One possibility to explain this would be 
that DLAs reach the primary plateau only for [0/H] > — 1. 

In the lower panel of Fig. [71 we plot (N/0) vs 
(0/H) -I- 12 for DLAs (squares) and local measurements 
(dots and crosses). The latters are from dwarf- irregular 
galaxy HII regions (van Zee & Haynes 2006), HII regions 
in spiral galaxies (van Zee et al. 1998), metal-poor emis- 
sion line galaxies (Nava et al. 2006; Izotov et al. 2006). The 
same remark as above is of course to be made here. It is 
apparent that DLA measurements are all below the local 
measurements. If a plateau is to be seen in DLA measure- 
ments, this is at (N/0) ^ —1.75. Other points are pretty 
well scattered in the diagram and we find difficult to con- 
firm the presence of a second plateau in Fig. [71 as claimed 
by Centurion et al. (2003). 

6. Conclusion 

We have measured the [0/Fe] and [N/0] abundance ratios 
in 13 DLA systems in the range —2.7 < [Fe/H] < —1 and 
added three systems from the literature. 
We show that the scatter in the [0/Fe] measurements is 
small around a mean value of [0/Fe] ~ 0.32±0.10. This 
small scatter is probably evidence for similar nucleosyn- 
thetic histories and efficient mixing of the gas. In the same 
range of metallicity, [0/Fe] in metal-poor stars is found to 
be about 0.7 dex if the measurements are not corrected for 
3D effects and 0.4 dex if corrected (Cayrel et al. 2004). 
In the [N/0] versus [0/H] diagram, we find that about a 
third of the measurements are located close to but below the 
so-called local primary plateau. If a plateau is to be defined 
for DLAs, its position is at [N/0] ^ —0.9 that is lower than 
the plateau measured locally in metal-poor dwarf galax- 
ies, [N/0] ~ —0.57 to —0.74. All our measurements are 
located below the local plateau which is not the case for 
measurements in metal-poor halo stars (Spite et al. 2005): 
about half of the stars have [N/0] > —0.9. However this 
may be due to small number statistics as a quarter of the 
DLA systems compiled by Henry & Prochaska (2007) have 
a [N/Si] ratio larger than the value corresponding to the pri- 
mary plateau for these elements. We find difficult to confirm 
the presence of a second plateau at lower value as claimed 
by Centurion et al. (2003) although the [N/0] distribu- 
tion is probably double peaked around [N/0] ~ —0.9 and 

[N/0] 1.35. This bimodality in the [N/0] distribution 

could be due to small number statistics as well although a 
possible bimodality in the [N/a] distribution has already 
been noticed by Prochaska et al. (2002) and Centurion et 
al. (2003); see also Henry & Prochaska (2007). 
Oxygen is mainly produced in short-lived massive stars 
and released into the ISM by Type II Supernovae explo- 
sions. Nitrogen is mainly produced in long-lived interme- 
diate mass stars and ejected in the ISM by stellar winds. 
This implies delays in the ejection of Nitrogen into the ISM 
even if it is primary. This is probably why most of the DLA 
measurements are below the local primary plateau. We note 
that the approximately constant observed [O /Fe] ratio does 
not support the suggestion by Pettini et al. (2002) that Iron 
should have the same evolutionary time-scale as Nitrogen 
so that DLAs deficient in Nitrogen should also be deficient 
in Iron. However, their argument holds only for part of the 
Iron production, the part produced by Type la Supernovae 
that have the same time-scale as intermediate-mass stars. 
Indeed, models with constant star-formation over a large 
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Fig. 7. Upper panel: Nitrogen to Oxygen abun- 
dance ratio, [N/0], versus Oxygen abundance [0/H]. 
[X/H] = (X/H)-(X/H)0. The range for the local primary 
plateau is indicated as two horizontal dashed lines at 
[N/0] = -0.57 and -0.74. The other dashed hue indicates 
the expected position in case of secondary production: the 
line is the extrapolation at low [N/0] values of local mea- 
surements. Lower panel: Same as above using the absolute 
abundances (X/H) = log N{X) - log iV(H) with X = O 
or N. The dots and crosses indicate local measurements 
from dwarf-irregular galaxy HH regions (crosses: van Zee 
& Haynes 2006), HH regions in spiral galaxies (crosses: van 
Zee et al. 1998), metal-poor emission line galaxies (dots, 
Nava et al. 2006), HH regions in Blue Compact Dwarf 
galaxies (Izotov & Thuan 2004) and emission-line galaxies 
from SDSS (dots, Izotov et al. 2006). 



range of duration produce a [Si/Fe] ^ 0.3 ratio about con- 
stant over a large metallicity range (Henry & Prochaska 
2007). These models cannot reproduce the scatter in the 
measurements but one could easily claim that constant star- 
formation rate is a simplistic assumption and that scatter 
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Fig. 8. N I and O i absorption lines in the DLA sys- 
tems: Q 0102-190, z = 2.926; Q 0112-306, z = 2.418; 
Q 0841+129, z = 2.476 and Q 0913+072, z = 2.618 from 
top to bottom respectively. Model fits are overplotted. 



p. Petitjean et al.: Nitrogen and Oxygen Abundance Ratio in DLA Systems 



9 





1 

0.5 






- nJ XI 134 


9 - 


nn^ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 i_ 


= I 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 





1 

0.5 



J 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1. 


- Si II A 13a 




1 I. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




-200 -100 100 200 -200 -100 100 200 
Relative velocity (km s~') 




-200 -100 100 200 -200 -100 100 200 
Relative velocity (km s~') 




-200 -100 100 200 -200 -100 100 200 
Relative velocity (km s~') 




I I I I I I I I I I I I I I I I I 



N I A 1199 

I I I 11 I I II i I 11 I I I I I I 



, N I A 1200:2 
hi i r n 11 I i T i I I I I I I I 





U I I III I I I M I I I I. 




I A 9501 

mi [mi 1 1 1 1 1 1 1 1 1 1 



-200 -100 100 200 -200 -100 100 20C 
Relative velocity (km s~') 

Fig. 9. N I and O i absorption lines in the DLA sys- 
tems: Q 1108-077, z = 3.608; Q 1337-^113, z = 2.508; 
Q 1337-hll3, z = 2.796 and Q 1340-136, z = 3.118 from 
top to bottom respectively. Model fits are overplotted. 
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Fig. 10. N I and O i absorption lines in the DLA sys- 
tems: Q 1409-^095, z = 2.456; Q 1409+095, z ^ 2.668; 
Q 2059-360, z = 3.083 and Q 2332-094, z = 3.057 from 
top to bottom respectively. Model fits are overplotted. 
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can arise from different star-formation histories (see also 
MoUa et al. 2006). 

As a consequence, considering an isolated galaxy, the 
[N/0] ratio could decrease sharply during a star burst when 
Oxygen is released by SNe (e.g. Contini et al. 2002). The 
scatter in the [N/0] ratio measurements could then be due 
to the intensity of the different bursts of star formation, ex- 
plaining the different [N/0] values. After some delay, cor- 
responding to the lifetime of intermediate mass stars, the 
[N/0] ratio could increase because of the release of most 
of the Nitrogen. The bimodality of the [N/0] distribution 
and the larger number of systems with low [N/0] ratio are 
a consequence of the delay between the releases of Oxygen 
and Nitrogen, because the duration of the release is small 
compared to the life-time of the stars (see also Ledoux et 
al. 2006b) . It is clear that increasing the sample size would 
be important to improve our knowledge on these issues. 
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